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1. Introduction 

Photonic bandgap (PBG) fibers guiding light in a hollow core surrounded by a 
cladding structure with a bandgap at a refractive index below the light line have at- 
tracted considerable attention since their first experimental demonstration by Cregan 



in the silica regions of the fibers. These findings suggest that dispersion 



et alP Such fibers have been proposed as candidates for highly linear and possibly 
low-loss transmission fibersP devices for particle transport P dispersion compensa- 
tiorP and gas nonlinearity experiments.^ In contrast to conventional fibers, the useful 
wavelength range is not limited by the absorption loss and nonlinearity of the base 
material. The recent fabrication of silica-based airguiding PBG fibers with attenua- 
tion coefficients below 30 dB/km over a considerable wavelength rang c 5 opens up for 
a wide range of practical applications. 

It is usually assumed that the influence of the base material on elementary fiber 
properties are negligible for this class of fibers, however the theoretical investigations 
performed up to now^' have primarily focused on establishing the shape and trans- 
mission windows of the guided modes, and have not, except in the case of circular 
Bragg fibers,^ provided a detailed modeling of key quantities such as group velocity 
dispersion (GVD), nonlinear coefficients etc. The purpose of the present work is to 
model two instances of a simple and well established design of airguiding PBG fibers 
with particular emphasis on the interaction between light and base material, which is 
here assumed to be silica. Specifically, we derive the fraction of the field energy present 
in the silica regions of the fiber, the nonlinearity coefficient (expressed as an effective 
area) arising from the material nonlinearity of silica, and the GVD including material 
dispersion effects. We demonstrate that the fraction of the field energy present in the 
silica is below 10% for both structures studied, and that the nonlinearity coefficients 
arising from silica are 1-2 orders of magnitude lower than what is obtained in the best 
silica-based large mode-area fibers. Furthermore, it will be shown that the GVD is 
considerably influenced by the dispersion of the base material, and demonstrate the 
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reason for this interesting effect by a detailed analysis of the material contributions 
to the GVD. 

The rest of the paper is organized as follows: In section 2, we describe the fiber 
designs to be investigated and briefly outline the theoretical approach adopted here, 
including the basic formulae for group velocity and nonlinearity coefficients. In section 
3, our numerical results are presented and discussed while section 4 summarizes our 
conclusions. 

2. Theoretical approach 

The two fiber structures to be investigated are both based on a cladding structure 
consisting of a triangular lattice of airholes, with a core defined by a larger airhole. 
The structure is characterized by three parameters: The distance between cladding 
hole centers, A, which is commonly denoted the pitch, and the diameters of core and 
cladding holes. We fix the core hole diameter at 3d, where d is the cladding hole diam- 
eter, and investigate two designs with o?/A=0.88 and (i/A=0.95 respectively. A design 
similar to the former has earlier been modeled by Broeng and co-workers^ whereas 
the latter resembles a low-loss airguiding fiber recently fabricated by Venkataraman 
and co-workers.^ A schematic picture of the core and nearest cladding region in the 
design with d/A=0.88 is shown in Fig. ^ 

In the present work we solve Maxwells equations by expanding the dielectric 
function and magnetic field vector in plane waves using a freely available software 
packageP Having obtained the magnetic field vector the electric fields are straight- 
forwardly calculated by use of Amperes law. The adoption of a planewave basis ne- 
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cessiates the use of periodic boundary conditions, however the interaction between 
nearest-neighbor repeated images of the guiding defect can be minimized by a proper 
choice of the transverse Bloch wave vector. ^ We use a supercell consisting of 8x8 ele- 
mentary cells of the triangular lattice comprising the cladding. The Fourier grid used 
for the plane- wave expansion has 64x64 meshpoints in each elementary cell for the 
structure with <i/A=0.88, and 96x96 meshpoints for the structure with <i/A=0.95. 
With these parameters, the dispersion coefficients and all other results, are converged 
within a few percent. 

The nonlinear coefficient of a fiber expresses the change in effective index of the 
guided mode arising from nonlinear effects for a given input power. The dependency 
of the nonlinear coefficient on the form of the guided mode is usually expressed by 
an effective areaP^ 
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Here P is the power launched into the fiber and nf is a material nonlinear co- 
efficient (related to the third-order nonlinear susceptibility) in units of W/m 2 . For 
conventional, all-silica, fibers A e g may be expressed asP^ 
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We have recently shown that for situations in which a substantial part of the field 
propagates in air the above definition must be generalized toP^ 



-4eff 




SiOa 
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Note that the integration in the denominator is now restricted to the silica parts 
of the fiber. This formula has been derived without making assumptions about the 
field energy distribution and is therefore applicable even in the extreme case of air- 
guiding PBG fibers. Of course, the A e g values obtained for these fibers have little 
to do with the physical extent of the guided modes, however the expression of the 
nonlinear coefficient in this form facilitates the comparison with more conventional 
fibers guiding light in silica or other materials. 

The GVD coefficient, D, is defined as: 
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where v g is the group velocity: 
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In the present case, where the dielectric function is piecewise constant, the group 
velocity in the presence of material dispersion effects, may be writtenP 
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where Ed is the fraction of the electric-field energy present in the dielectric and 
Vg is the group velocity in the absence of material dispersion. The latter may be 
calculated directly from the fields as!^ 

_ Re([E* x HP 
9 (H,H) 1 ' 



Thus, the group velocity v g can be evaluated directly from the fields once the 
guided mode has been obtained, and the dispersion coefficient can then be calcu- 
lated by a numerical first-order derivative. This procedure requires that u>, Ed and 
Vg are evaluated at the silica refractive index appropriate for u, which in the present 
work is achieved by a self-consistency procedure.^ The self-consistent calculations are 
compared with calculations assuming a fixed value of the silica refractive index, n, 
in order to assess the importance of material dispersion effects. In the selfconsistent 
calculations we use the Sellmeier formula for the frequency dependence of the silica 
refractive index, with the coefficients reported by Okamoto.^ 
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3. Numerical results 

In this work we focus on the guidance of the fundamental mode (whose major trans- 
verse part is circularly symmetric) in the lowest bandgap. Initially, we will consider 
the case of a fixed silica refractive index n=lA5. For this value of n, the fiber with 
g?/A=0.88 is found to have a narrow transmission window for the fundamental mode 
between A/A=0.724 and A/A=0.685, whereas the fiber with c//A=0.95 has a some- 
what wider transmission window between A/A=0.617 and A/A=0.533. As will become 
clear later these transmission windows show some dependence on the material refrac- 
tive index, which translates into a dependence on the physical value of the pitch 
(since this controls the physical wavelength of the light in the guided mode). In both 
fiber designs higher-order modes are present in part of the transmission range of the 
fundamental mode. For g?/A=0.88 we find that second-order modes are present in the 
fundamental bandgap in the lower three-quarters of the transmission window for the 
fundamental mode. For d/A=0.95 the second-order modes leave the bandgap some- 
what earlier, when the fundamental mode is roughly in the middle of the bandgap. 
Since the question of determining the single-mode wavelength regions of the fibers is 
complicated by the possibility of guidance in the higher-order bandgaps, and is not a 
primary concern in this paper, we have not attempted a precise determination of the 
transmission windows for the second-order modes. 

In Fig.[2fa) the fraction of the electric field energy present in the silica part of the 
fibers (Ed in Eq. (jUJ)) is plotted as a function of the distance between the frequency 
of the fundamental mode and the lower band-gap edge normalized to the gapwidth. 
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Both fiber designs show the same qualitative behaviour: Ed rises as the mode enters 
or leaves the gap, and therefore a minimum is present inside the transmission window. 
However, for the design with d/A=0.88 the minimum is present in the low-frequency 
part of the transmission window, whereas for c?/A=0.95 the minimum is shifted close 
to the high-frequency transmission edge. It is also noteworthy that the frequency 
derivative of E^ is quite large, since the transmission windows are narrow. This has 
important consequences for the dispersion properties of the fibers. 

In Fig. |2Jb) the effective areas, as calculated from Eq. J3J), are plotted for the 
two fiber designs. The results for d/A=0.88 have been multiplied by a factor of 10 to 
facilitate comparison. As expected, very large A e g values are found, signifying very low 
nonlinear coefficients. In index guiding microstructured fibers in the large-mode area 
regime (A <C A) the effective area A e g ~ a x (d/A)~ 1 A 2 with a numerical prefactor 
a of the order 0.5.^ The fibers are typically operated close to the endlessly-single 
mode limit (d/A ~ 0.45) so that A e g ~ A 2 . Typical values of A are 10-20 /im, so 
that A e g ~100A 2 for A ~l/im. Thus, the present results for airguiding PBG fibers 
indicate a lowering by 1-2 orders of magnitude of the nonlinear coefficient compared 
to typical index-guiding large-mode area microstructured fibers available. Still, it is 
interesting to observe the significant variation of A eS over the transmission window, 
and the strong dependence on cladding design of the nonlinear coefficients. In the 
fiber with d/A=0.88 a decrease of the effective area with increasing frequency is seen, 
corresponding to the increasing fraction of field energy in silica (see Eq. (j3J)). For 
d/A=0.95 the opposite trend occurs, due to the shift of the minimum in Ed. Of 
course, the effective areas reported here relate to the nonlinearity coefficients and 
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have little to do with the physical size of the modes. This is better estimated from 
the standard definition of effective area, Eq. (J2J), which for both fiber designs is found 
to be comparable to the area of the hollow core, indicating that the guided mode is 
well localized. 

In Fig. El GVD results for three different fiber designs are reported. For the 
fiber with d/A=0.88 we have investigated two values of the pitch, A=0.8 /im and 
A=2.4 /im. For the design with <i/A=0.95 we show results for A=1.0 /im. Both the 
results of self-consistent calculations and of calculations with a fixed value of the silica 
dielectric constant are shown. It can be seen that a change in the silica refractive 
index, n, shifts the transmission windows, and thereby the dispersion curves. Due 
to the steepness of the dispersion curves this implies that the dispersion at a given 
wavelength is strongly dependent on n. Therefore, the waveguide GVD calculated 
at n=1.45 (the solid curves), which is the refractive index of silica at a wavelength 
of 1.05 fim, gives a poor prediction of the true chromatic dispersion (as given by 
the self-consistent calculations, reported by the dotted curves) at other wavelengths. 
The agreement is considerably improved by choosing a fixed index suitable for the 
wavelength of the guided mode. The results of such calculations are reported by the 
dashed curves. However, there is still a noticeable difference between the dispersion 
curves calculated at a fixed n and the self-consistent results. The differences are of 
the same order of magnitude as the material dispersion of homogeneous silica at the 
wavelengths in question, and are seen to change sign over the transmission window. 
In Fig. Of a) it is interesting to notice that the dashed and dotted curves do not tend 
towards each other at the shortest wavelengths of guidance even though the value of 
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71=1.46 used for calculation of the dashed curve corresponds to a wavelength of ~550 
nm for pure silica. Instead, the curves cross at somewhat longer wavelength. These 
findings indicate that material dispersion effects play a significant role despite the 
small percentage of field energy present in silica. 

In order to obtain a more detailed understanding of the influence of material 
dispersion, we return to Equations and (0). The derivative of the group velocity 
with respect to frequency may be written: 



dv g _v g ( 1 dVg dVg de /E d d\ne oj fdE d d\iae cPlne" 
~du~ = U£~d/3 + ~de~d^~ Vg \T~cLj~ + 2 I ~du7~du7 + d du 2 



In this formula, d/du (d/d/3) denotes a derivative with respect to uj (j3) for a 
fixed value of e, whereas d/du (d/d/3) denotes a derivative including the variation of 
e with uj (and thereby (3). If dispersion in the base material of the fiber is neglected 
only the first term contributes. Using the equationP 



dv° g E d u 8E d 



de 2e 9 2e 8(3 



and approximating dE d /duo m dE d /du, which we have found to be reasonably 
well justified even for the airguiding fibers studied here, we can write: 
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Thus the GVD may be separated into a part independent of material dispersion 
effects (first term in Eq. (|lUp). a part proportional to Ed and a part proportional to 
the frequency derivative of E d . Using Eq. (£Q) the GVD is found to be: 



where is the GVD in the absence of material dispersion, but evaluated at the 
silica refractive index appropriate for the uj value in question. In Fig. 0] the material 
dispersion D mat defined by Eq. (jTTj) is plotted for the two designs guiding at short 
wavelengths. For the design with <i/A=0.88, where 7-9% of the field energy is in silica, 
the material dispersion ranges between ~0 and -200 ps/nm/km, whereas for the design 
with d/A=0.95, and only 2-3% of the field energy in silica, the material contribution 
to the GVD ranges between -50 and 50 ps/nm/km. Also shown in Fig. 0] is the 
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difference between the waveguide dispersion at fixed, suitably chosen, n (dashed lines 
in Fig. EJ) and the self-consistent dispersion coefficients (dotted lines in Fig. EJ). Exact 
correspondence between solid and dashed curves in Fig. |U is not to be expected since 
the waveguide dispersion, D^ c in Eq. (JTTjl is evaluated at the self-consistent value of 
n, however, it can be seen that the major part of the discrepancy between Dw and 
D sc in Fig. |3] can be attributed to the intrinsic material dispersion as expressed by 
Dmat- Since the material dispersion of homogeneous silica in this wavelength range is 
between -250 and -400 ps/nm/km, the D mat values reported in Fig. H|are surprisingly 
large considering the small values of Ed- 
it is evident from Eq. (JTTJ) that the contribution of material effects to the total 
GVD of a fiber is composed of a part proportional to Ed and a part proportional to 
d Ed /duo. Herein lies the origin of the surprisingly large D mat values for the airguiding 
fibers: Although E d is small in these fibers, as is evident from Fig. 12(a), this is not 
the case for dEd/duj. In Fig. EJwe plot the ratio, Rb between the third and second 
term in Eq. ([lip: 
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The results for the airguiding fibers in Fig. OJa) are compared to results for index 
guiding microstructured fibers, reported in Fig. El[b). The latter have a triangular 
cladding structure similar to the airguiding fibers, but with a solid silica core defined 
by a missing airhole. It can be seen that | Rd | for the index guiding fibers is every- 
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where below unity, even in the rather extreme case of d/A=0.8, A=0.34 fim, where 
~15% of the field energy is located in the airholes. For the airguiding PBG fibers 
| Rd I is 1-2 orders of magnitude larger. Thus, for index-guiding microstructured 
fibers the main contribution to material dispersion effects comes from the second 
term in Eq. (jlljl . whereas for airguiding PBG fibers the contribution from the third 
term dominates. 

Because airguiding PBG fibers have the major part of the field energy propagating 
in air, the choice of base material is less limited by requirements of low loss and/or 
nonlinearity than is the case for standard fibers, or index-guiding microstructured 
fibers. On the other hand, the results presented in this work show that the dispersion 
properties of the material may still have a significant impact on the total GVD of the 
fiber. These observations suggest that dispersion engineering through the choice of 
base material may be a possibility in these fibers. A simple example of the possibilities 
is shown in Fig. [HI The usual three-term Sellmeier polynomial describing the material 
dispersion of silica has been modified by a fourth term describing the addition of an 
(dopant) absorption line close to the transmission window of the fiber. The modified 
Sellmeier polynomial reads: 
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with A4 chosen as either 0.5 /im (P500 in Fig. |HJ) or 0.640 /im (P640 in Fig. EJ). 
It is evident that significant shifts of the dispersion curve can be obtained simply 
by addition of (impurity) absorption centers to the silica matrix. A more general 
approach would of course be to vary the composition of the base material itself as 
could readily be done in, e.g., polymer fibers. Such dispersion engineering could, 
for instance, be of interest for fibers applied to the kind of gas-phase nonlinearity 
experiments whose feasibility was recently demonstrated by Benabid and co-workersP 

4. Conclusions 

In conclusion, we have investigated various aspects of the interplay between base ma- 
terial and the fundamental guided mode in silica-based airguiding PBG fibers. For 
the two designs studied here, between 2 and 9% of the electric field energy was found 
to reside in the silica parts of the fiber. The nonlinearity coefficient was expressed in 
terms of a generalized effective area, which was found to be 1-2 orders of magnitude 
larger than what can be obtained in index-guiding microstructured fibers. The influ- 
ence of material dispersion on the total GVD of the fibers was investigated and was 
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found to be of the same order of magnitude as in other fiber types having most of the 
field energy residing in silica. This effect was traced to the fact that the variation with 
frequency of the field energy in silica is much more rapid in airguiding PBG fibers 
than in other fiber types. These results suggest that dispersion engineering through 
the choice of base material may be an interesting possibility in airguiding PBG fibers. 
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List of Figure Captions 

Fig. 1. Schematic picture of one of the structures (with d/A=0.88) under study. The 
black circles are airholes, while the white areas are the silica regions. Only the core 
and innermost cladding region is shown. 

Fig. 2. Field energy fraction in silica (a) and effective calculated from Eq. 

(JHJ) (b). The effective area curve for the fiber with d/A=0.88 (solid curve in (b)) has 
been multiplied by 10 to facilitate comparison. 

Fig. 3. Dispersion curves for the fundamental guided mode of three airguiding 
PBG fibers with various values of rfand A. (a): d/A=0.88, A=0.8 /im. (b): d/A=0.88, 
A=2.4 /im. (c): d/A—0.95, A=1.0 /im. Solid curves report waveguide dispersion (Dw) 
calculated at n=1.45, dashed curves report waveguide dispersion at values of n suit- 
able for the wavelength interval spanned by the transmission window, and dotted 
curves denote results of self-consistent calculations. 

Fig. 4. Material dispersion D mat (solid lines), defined in Eq. (JHJ), for the two fiber 
designs guiding at short wavelengths. The dashed lines report the difference between 
the dotted and dashed curves in Fig. 3 for comparison. 

Fig. 5. Plots of the quantity Rd, defined in Eq. (JHJ), for two airguiding PBG 
fiber designs (a) and two index-guiding fibers (b) having a cladding structure similar 
(although with smaller airholes) to the airguiding PBG fibers. 

Fig. 6. Dispersion curves for two fibers with added absorption resonances in the 
base material at either A=500 nm (P500) or A=640 nm (P640) compared to the 
undoped result. 
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